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Abstract. We demonstrate universal quantum control over chains of ions in 
a surface-electrode ion trap, including all the fundamental operations necessary 
to perform algorithms in a one-dimensional, nearest-neighbor quantum computing 
architecture. We realize both single-qubit operations and nearest-neighbor entangling 
gates with Raman laser beams, and we interleave the two gate types. We report average 
single-qubit gate fidelities as high as 0.970(1) for two-, three-, and four-ion chains, 
characterized with randomized benchmarking. We generate Bell states between the 
nearest-neighbor pairs of a three-ion chain, with fidelity up to 0.84(2). We combine 
one- and two-qubit gates to perform quantum process tomography of a CNOT gate in 
a two-ion chain, and we report an overall fidelity of 0.76(3). 


RAGS numbers: 03.67.Lx, 03.67.Bg, 32.80.qk 


quantum computing, quantum control, trapped ions, surface-electrode trap 

1 . Introduction 

The size of the Hilbert space describing a quantum computer grows exponentially with 
the number of qubits. Despite this, any quantum computation can be implemented to 
arbitrary accuracy by composing elements from a set of one- and two-qubit operations 
DEI- Any such operation set is called a universal gate set. Universal gate sets have 
been demonstrated with NMR [3] , superconducting mu, neutral atom [6] , and trapped 
ion mu qubits. 

Trapped ions are attractive for building an extensible quantum system due to their 
long coherence times and ability to be entangled via photons and phonons PE]. Ion 
trap quantum computer architectures could use transport [n] or photonic interconnects 
[12] to provide connectivity between disparate regions of the quantum computer. These 
architectures can benefit from the fabrication precision and scalability of microfabricated 
surface-electrode ion traps [13]. While components of a universal gate set have been 
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shown in microfabricated traps (see e.g. [a OEl dSl El [13), a complete nniversal gate 
set has not yet been demonstrated in a single system. 

In this work, we demonstrate a nniversal gate set on two- and three-ion chains in a 
snrface-electrode trap. Our architecture is based on transporting the chain sequentially 
through a single gate location to provide ion addressing; the gate beams are pulsed on 
only when the chain is stationary. We perform both single-qubit operations and nearest- 
neighbor entangling gates. We measure single-qubit gate hdelities for chains of two to 
four ions using randomized benchmarking, Bell state hdelities in the nearest-neighbor 
pairs of a three-ion chain, and process hdelity of a cnot gate in a two-ion chain. The 
process hdelity measurement demonstrates our ability to interleave one- and two-qubit 
operations in the same experiment. All features are scalable to longer ion chains. 


2. System overview 


The experimental system is shown schematically in hgureiHa). We trap chains of 
ions 60 pm above the surface-electrode ion trap described in [T9]. This trap features 
through-chip vias instead of wirebonds. The improved optical access provided by this 
feature allows for tight focusing of the Raman gate beams, which propagate parallel 
to the surface of the trap. Qubits are encoded in the hyperhne clock states of 
with |0) = ^51/2 \F = 0,mF = 0) and |1) = |F = l,mF = 0). Details on trapping, state 
preparation, and readout of ^^^Yb’*' ions can be found in [20] . 



Figure 1. (a) Schematic diagram (not to scale) of control and individual detection. 
Both single and two-qubit gates are performed with the same Raman gate beams, 
which address a pair of ions. The quantization field B is parallel to the Raman beams, 
(b) With fixed Raman beams, the entire chain is transported along the trap axis to 
address different ion pairs in a “cascaded” fashion. We apply the series of gates Ut to 
generate an arbitrary unitary Ui on each ion i. 

We load individual ions and merge each ion into a harmonic well. In this 
way, we deterministically build chains of a desired number of ions. During state 
































Universal Control of Ion Qubits in a Scalable Microfabricated Planar Trap 


3 


detection, each ion in the chain is imaged onto a separate channel of a multi-channel 
photomultiplier tube (PMT) |2l]. This allows us to detect the state of each ion in 
the chain independently. Despite efforts to reduce crosstalk by using only every fourth 
PMT channel, we observe 5% crosstalk between ions due to coma in the imaging system. 
This crosstalk is calibrated and accounted for with a joint determination of which ions 
were bright in a given experiment set. With 400 ps of fluorescence collection, our state 
detection hdelity is 0.98(1). 

A single pair of Raman beams provides both one- and two-qubit gates. The beams 
are derived from a mode-locked tripled YAG laser, in a conhguration similar to that 
described in [221 [23]. The circularly polarized beams propagate perpendicular to the 
trap axis and parallel to the quantizing magnetic held. In this conhguration, the gate 
beams do not couple to the axial motional modes or the \F = l,mF = ±1) Zeeman 
levels. Furthermore, the trapping potentials are designed so that one set of radial secular 
modes is rotated normal to the beam direction. This limits coupling of the Raman beams 
to only a single set of motional sideband transitions along with the carrier, reducing the 
ehect of oh-resonant shifts and removing the dependence of the carrier transition on 
the temperature of these modes. After Doppler cooling all modes, we perform sideband 
cooling with the Raman beams only on the radial mode used for gates, resulting in 
sub-quanta occupation of excited motional states. 

Mplmer-Sprensen (MS) gates between pairs of neighboring ions provide the 
entangling gates for this system. The Raman beam in one direction (consisting of the 
red- and blue-detuned beam^ is focused to a waist (1/e^ half-width of the intensity 
prohle) of roughly 7 pm while the opposing Raman beam has a waist of 13.5 pm. The 
average power in these beams is on the order of 10 and 100 mW, respectively. The 1/e 
half-width of the Raman interaction Rabi rate (see hgure [2|) is approximately equal to 
our ion spacing of about 6 pm, permitting pairwise addressing of ions within a chain as 
in hgure m The Raman beams remain static throughout the course of the experiment; 
we transport the ion chain to bring different pairs of ions into the beam. 

Still addressing ions pairwise, we employ the cascading transport scheme shown in 
hgure [T](b) to implement an arbitrary set of one-qubit unitary operations {f/*}, where i is 
the ion index. The ion chain is hrst transported such that the i = 0,1 ion pair (i.e. the 01 
pair) is illuminated by the laser. The unitary Uq = Uq is performed, leaving the endmost 
ion in the desired state. Then the ions are transported to address the subsequent pair. 
The operation Ui = UiUl results in the net operation UiLIq = UiUqUq = U\ on the 
second ion, accounting for the ehect of the hrst operation. This procedure is repeated 
sequentially along the chain until all ions are targeted, resulting in the desired gate on 
each ion. Universal control is achieved in this transport architecture by interleaving 
such cascaded single-qubit gates with pairwise MS gates. 

In this work, ions are only transported over very short distances corresponding to 
the spacing between ions or the length of the chain. Each transport step is performed 

f These beams are generated with a single acousto-optic modulator, the output of which is imaged 
onto the trapping region with several telescopes to ensure the two beams overlap at their foci. 


Universal Control of Ion Qubits in a Scalable Microfabricated Planar Trap 


4 


adiabatically in 100 /is. We varied the transport dnration between 100 ps and 1 ms bnt 
could not see any effect on gate fidelity at the present error levels due to heating. In 
addition to shuttling the ions through the gate beams, we adjust the strength of the 
axial harmonic well to optimize the ion spacing for each type of gate. For two-qubit 
operations, where a strong coupling between the ions’ motion is desired, the axial well 
is stronger than for single qubit operations, where crosstalk to unaddressed ions needs 
to be minimized. For Doppler and sideband cooling, the typical axial trap frequency 
is 271 X 0.5 MHz, and the radial frequencies are 27 r x (1.8, 2.1) MHz. The lower radial 
frequency corresponds to the mode addressed for all gates. We reduce the axial frequency 
by a factor of y/2 from the cooling trap frequencies to spread the ions out for single-qubit 
gates, and we increase it by \/L5 times to further split the radial modes for two-qubit 
operations. Similarly, for detection the ion spacing is optimized to match the spacing 
of the photomultiplier array. 
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Figure 2. Characterization of the intensity profile and phase front of the Raman 
control beams. The position of each ion in a 4-ion chain is shown across the top 
when ions 1 and 2 are addressed. The beam power was calibrated for tt/2 rotations 
on the addressed pair. Open (filled) circles use bare (PBl stabilized) pulses, and 
markers are connected for clarity. The PBl pulse sequence stabilizes Rabi rate errors, 
resulting in the flatter slope at rotation angle 7r/2 around ions 1 and 2. Additionally, 
for bare rotation angles less than 7r/4, the PBl sequence suppresses rotations. The 
slight asymmetry between ion 0 and 3 is due to beam coma. 


3. Single-Qubit Control 

As described above, all single-qubit control is implemented via pairwise addressing. For 
each ion pair, we calibrate the ion chain position to balance the carrier transition Rabi 
rate on the two addressed ions. We then calibrate the gate beam intensity to give a 7r/2 
rotation in 4 ps. These calibrations are stabilized against Rabi rate errors using PBl 
compensating pulse sequences [2ll [25l [26], which are composed of the calibrated 7r/2 
rotations. Including hardware programming overhead, a PBl composite 7i/2 operation 
is performed in 120 /rs. Arbitrary single-qubit unitaries are constructed from these 
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composite 7r/2 rotations about different axes in the equatorial plane of the Bloch sphere. 
The angle the axis of rotation makes with the x-axis is set by the phase of the Raman 
beams and is controlled via an acousto-optic modulator. We developed a compiler (see 
to efficiently implement arbitrary unitaries in the fewest 7r/2 rotations. 

Implementation of an arbitrary set of unitaries across the chain under actual 
experimental conditions is more complicated than the idealized scheme presented in 
£gure[T](b). First, wavefront curvature across the gate beams causes ions to experience 
a position-dependent phase of the Raman interaction (see figure Ej), resulting in a 
spatially varying 7r/2 rotation axis. Second, the pairwise addressing of our gate 
beams is imperfect, and ions adjacent to the addressed pair are observed to flop at 
reduced but signihcant rates. To a large degree, this crosstalk is mitigated by our 
choice of compensating pulse sequence since PBl, a so-called “passband” sequence, 
also suppresses small excitations (figure [2]). Due to comatic aberrations, the crosstalk 
is asymmetric, and it is a reasonable approximation to treat the right-hand side of 
our gate beam as crosstalk-free when PBl is used. Therefore, we only calibrate and 
correct for the small rotation angle of left-side crosstalk. This asymmetry dictates the 
direction in which we cascade single-qubit operations. As an example of independent 
single-qubit control, we prepare each ion in a four ion chain in all 3*^ combinations of 
{|0).|l).l+) = (|0) + |l))/^/2} (agureO, 


Appendix A) 



Figure 3. As a demonstration of independent single-qubit control through pairwise 
addressing, we perform all combinations of the set {I,H,X} on each qubit in a 4-ion 
chain after initializing to |0000). Probabilities are offset one unit for each ion. Each 
unique preparation is labeled by an experiment index and was repeated five times to 
show the spread in results. 

To characterize the single-qubit gate hdelity, we employ randomized benchmarking 
|27] . Starting with all ions in state |0), a randomly chosen Clifford gate is performed 
on each ion, cascading the gates across the chain as described above. Each Clifford 
operation requires on average 400 ps of Raman pulses and is followed by 100 ps transport 
between gate locations. Subsequent passes build a random sequence of Clifford gates 
of length L on each ion. The hnal pass consists of an inverting gate calculated to take 
each ion to |1). The results for one ion in a two-ion chain are shown in figure 111(a) for 
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varying sequence lengths. We performed benchmarking on two- to four-ion chains and 
£t the average gate fidelity Frb for each ion, following [28], to 

f'RB = i + i(l-2e„)(l-2e,y. (1) 

where Eg is the average gate error and Em includes error from state preparation and 
measurement as well as error in the inverting gate. The average Em for hts on all 
chain lengths is 0.025(2). We also performed a three parameter £t as in [29], where the 
constant offset is allowed to vary from I, however these hts gave non-physical results 
with negative Em- Because the measurement imbalance in our apparatus is small, the 
two-parameter fit with offset hxed at | is justihed; a typical reduced-y^ for the hts is 
1 . 6 . 




Figure 4. (a) Typical randomized benchmarking data, here for ion 1 in a 2-ion chain. 
Gray crosses show the data spread at each sequence length, and the blue circles are 
the average data with error bars reflecting the standard error in the mean. The fit is 
described in the text, (b) Average Clifford gate error for each ion as a function of ion 
chain length. Error bars represent the statistical error of the fit. 


As can be seen in figure |D(b), we measure average single-qubit gate fidelity as high 
as 0.970(1) for two-ion chains. However, the average gate error increases with ion chain 
length and position within the chain in the cascade direction. Longer chains require 
more calibrations and more operations, which lead to more time for dephasing between 
successive operations on a given qubit. As a result, one expects gate error to increase 
with chain length. 
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That the end ions, i.e. ions 2 and 3, have a higher gate error can be explained by 
imperfect crosstalk correction. As we noted in the discussion of hgure [2l the wings of 
our comatic gate beams cause small-rotation crosstalk on ions adjacent to the addressed 
pair. The PBl sequence suppresses the average value of this crosstalk to the point where 
it is negligible on the left-hand (low ion number) side of the beam; on the right-hand 
side, we calibrate and correct for the remaining crosstalk in our cascaded gate sequences. 
Therefore, ions 0 and 1, unlike ions 2 and 3, are not subject to this crosstalk. Despite 
our careful calibration, we suspect that the increased slope due to PBl (hgure [2]) may be 
increasing our sensitivity to crosstalk huctuations on the right-hand (high ion number) 
side of the beam. Improving the passive stability of our gate beams should increase 
the overall single-qubit gate hdelity. We note that the average gate error for ion 2 
decreased from a three-ion to four-ion chain. This is most likely due to relatively poorer 
calibration; systematic errors are not included in the error bars. 

4. Two-Qubit Gates 

f.l. Two-Ion Entanglement 

To complete our universal gate set, we require an entangling gate. Here we employ the 
Mplmer-Sprensen (MS) gate on a pair of ions [30]. During MS gates, we strengthen our 
axial conhnement to increase the frequency splitting of our radial modes in comparison 
to single-qubit gates. This reduces motional excitation of undesired modes, yielding 
higher hdelity. State evolution of a typical gate is shown in hgure [5](a) where ion motion 
has been excited and de-excited twice per gate, i.e. two “loops” per gate, resulting in a 
169 /is gate time. This is substantially faster than our PBl compensated implementation 
of single-qubit unitaries. 

To verify entanglement following the MS gate. Bell state hdelity estimation is 
performed using the “parity hopping” technique [31]. The resulting parity curve is 
shown in hgure [5](b). Bell state hdelity is measured to be Feeii = 0.93(2). We expect 
this hdelity is largely limited by imperfect overlap of the red and blue sideband beams 
during the MS interaction; a redesigned optics system will remedy this issue. This MS 
gate, along with the individual qubit rotations in section (3] comprise our universal gate 
set for a two-ion chain. In section (5] we demonstrate combined one- and two-qubit gates 
to produce a cnot gate. 

4-2. Pairwise Entangling Gates on Three Ions 

Universal control of ion chains requires the ability to entangle any pair of ions. Within 
our architecture, we entangle neighboring ions in pairs. However, this entanglement can 
be transferred to non-neighboring pairs through SWAP gates, which can be composed 
from CNOT gates. In this section, we characterize entangling gates on pairs of 
neighboring ions within a three-ion chain, the last component of a universal gate set for 
three ions. 
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Figure 5. (a) Population evolution during the Mplmer-Sprensen interaction of the four 
two-ion states: Pqo, Pqi, Pw) Pii- (b) Parity measurement (Pqo + — Poi — Pw) of 

the entangled ions at gate time of 169 /xs. 


Table 1. Parity measurement for pairwise entanglement in a three-ion chain. Pqo and 
Pii are the measured populations for the state where the unaddressed ion is in |0), 
i.e. for an MS gate on the 12-pair, the states |000) and |011). Similarly, the parity 
amplitude is extracted with the unaddressed ion in state |0). Pi is a measure of the 
population build-up on the unaddressed ion, where we trace over the addressed ions. 


Ion pair 

-Poo 

Pii 

Parity amp. 

Pi 

PBell 

No echo 01 

0.32(1) 

0.45(1) 

0.66(1) 

0.18(1) 

0.72(1) 

Echo 01 

0.37(1) 

0.46(1) 

0.62(1) 

0.10(1) 

0.72(2) 

No echo 12 

0.47(1) 

0.41(1) 

0.81(1) 

0.094(7) 

0.84(2) 

Echo 12 

0.41(1) 

0.47(1) 

0.76(1) 

0.023(3) 

0.82(2) 


To estimate the MS gate fidelity, we measure F^ew as in the previous section. All 
ions are prepared in |0), the two targeted ions are shifted into the laser beam, then 
an MS gate and analysis operation are performed. Parity is calculated using only the 
state populations where the unaddressed ion remains in the desired |0) state, and the 
resulting fidelities are summarized for each targeted pair in table [H 

We observe moderate build-up of undesired population in the |1) state on the 
unaddressed ion during the MS gate; this is reported as Pi in table [H Due to the 
spatial asymmetry of our gate beam (shown in the rotation angle of figure [2]), the MS 
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crosstalk on the 01 pair is worse than the 12 pair. In Appendix B, we describe an 
“echoing” sequence which decouples a subset of ions from a general MS interaction. For 
an MS gate on two ions in a three-ion chain, we perform the MS gate with two motional 
loops. The echo is inserted after each loop and flips the spin of the untargeted ion (or 
equivalently the spins of the targeted ions) with a Y-gaXe. Early results show that the 
echo sequence reduces Pi; however, the Bell-state hdelity is not improved. This is likely 
due to the error introduced by the additional single-qubit operations used in the echo. 
We expect that increased single-qubit gate fidelity should make the echo more effective 
and allow us to perform high-hdelity MS gates on any adjacent pair of ions in longer 
chains. 


5. CNOT Characterization 


For a demonstration of a small circuit using our universal gate set, we perform the 
CNOT gate on two ions. The cnot is built from an MS gate wrapped with a particular 
set of single qubit rotations, as shown in £gure[ 6 ](a). In figure [ 6 ](b,c) we represent how 
the CNOT circuit fits into the cascading architecture discussed in section |2l Single-qubit 
rotations are performed in two cascading passes, one before the MS gate and one after. 
As in section 01 axial conhnement is increased for the MS gate and then subsequently 
relaxed for single-qubit operations, adding an additional 50 /rs for each. We utilize four 
steps of transport (100 /is each, as before), excluding the transport needed for state 
preparation and detection. 

Successful combination of one- and two-qubit operations requires accounting for the 
Raman phase difference due to increased axial conhnement during MS gates. We add z- 
rotations (P^) before and after the MS gate (with rotation angles 933 and ip^) to ensure 
it is in phase with the single-qubit rotations. In addition, we add the optimization 
operations shown in figure [ 6 ](a). These gates do not affect the outcome of the ideal 
circuit, but they allow us to minimize the number of physical operations in the compiled 
single-qubit cascades by varying pi and p 2 - Further discussion of this procedure is 


provided in Appendix A 


As a diagnostic tool, we measure the population evolution throughout the 
implementation of the cnot (markers in figure [ 6 ](d)). To measure the populations 
at a given point within the cnot sequence, we trigger an override on the gate beams, 
shutting down all subsequent operations. This freezes the populations, which are read 
out at the end of the cnot sequence. When combined with on-the-hy simulation of 
ideal rotations (points connected by solid lines), such a scan has proved a useful tool 
for diagnosing calibration errors. We have also built an optimization routine into our 
simulation tool to hud the best parameter calibrations to match observed populations. 
Thus, from a small set of carefully chosen gate scans like £gure[ 6 ](d), we can quickly 
improve overall fidelity. 

To systematically characterize our cnot gate, we perform quantum process 
tomography (QPT) [32]. We prepare and analyze each ion in all combinations of the set 
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Figure 6. (a) Circuit diagram to build a CNOT gate (dashed box, controlled by lower 
ion) from single-qubit rotations and an MS gate (orange box). After initialization to 
|00), the ions are prepared in |10) (yellow box), (b) Schematic of pairwise addressing 
for compiled CNOT. Each purple block represents a PBl pulse sequence stabilized 7r/2 
rotation, (c) Ion transport during CNOT implementation, (d) Measured population 
evolution throughout our compiled CNOT. Each numbered gate is a 7r/2 rotation 
except for the orange region in which the MS gate is shown in 2 steps. Simulated 
expectations are connected with solid lines and agree very well with the data. For 
input state |10), the CNOT gate yields the expected output state |11) with probability 
0.85. 

{|0), |1), 1+), |—)^ = (|0) — i |1)) /\/2} resulting in 256 measurements. The preparation 
and analysis gates are not compiled into the single-qubit gates in the CNOT itself. 
Rather, two additional passes of single-qubit operations are used, similar to the state 
preparation (yellow box) in hgure |6l Each measurement is repeated 160 times and 
the maximum likelihood probability to be in |00) is used to generate the measured 
process matrix Xmeas, shown in hgure [71 The hdelity of our operation with the ideal 
CNOT process, Xideai is measured as F = Tr(xideaiAmeas) = 0.76(3), where we have 
reported the mean and standard deviation of 10^ parametric bootstrap resamplings of 
our measurements. From the resampled Xmeas, we estimate a typical matrix element 
error bar to be 0.02. There are only 16 nonzero matrix elements of Xideah all in the real 
part. Nearly all other measured matrix elements are consistent with zero within three 
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Figure 7. Quantum process tomography of CNOT on a 2-ion chain. Bar height reflects 
magnitude of the real (imaginary) part in the left (right) plot with color giving the 
sign. Ideal matrix values are shown as colored wire frames. All imaginary elements 
should be identically zero. The apparent CNOT fidelity from QPT is 0.76(3). Figure 
data included as supplementary data on-line. 

Using a simple nncorrelated error model to acconnt for fonr single-qnbit gates (from 
randomized benchmarking fidelity) and one MS gate (from Bell state fidelity), we expect 
an approximate fidelity of {F-^BYFBeW = (0.97)^ x 0.93 = 0.82. Yet, single-qnbit errors 
in preparation and analysis also rednce QPT fidelity estimates. With the same QPT 
procednre, we measnre a two-qnbit identity operation hdelity of 0.95, leading ns to 
estimate a QPT hdelity of 0.78 in good agreement with onr measnrement. 

6. Outlook 

We have demonstrated a nniversal gate set for two- and three-ion chains for the hrst time 
in a snrface-electrode ion trap. Universality is realized throngh a novel control scheme 
which leverages the precision transport capability of onr trap. This scheme is extensible 
to longer chains. Importantly, we demonstrated nniversal control in a single experiment 
by implementing a circnit consisting of single-qnbit rotations and an MS gate to perform 
a CNOT gate. Single-qnbit gates on fonr-ion chains were characterized as well as echoed 
MS gates on pairs within a three-ion chain. Minor system improvements cnrrently 
nnderway shonld improve hdelities to levels where basic qnantnm algorithms are within 
reach. 

Onr resnlts bolster the promise of scalable ion trap qnantnm compnting. Throngh 
the integration of jnnctions, microfabricated ion traps have already been extended to 
two-dimensional architectnres. Integrated optics are a demonstrated path to a scalable 
photonic interface with ions in large trap arrays. We have shown that microfabricated 
traps are also compatible with high-hdelity nniversal qnantnm control of ion chains. 
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Appendix A. Compiler 

In our control scheme, we must implement arbitrary unitary operations on individual 
qubits for universal control. Even for algorithms designed around a discrete set of 
single-qubit gates, accounting for the crosstalk shown in hgure |2] would require unitary 
operations outside that set. Here, we describe how we construct arbitrary unitaries from 
our calibrated, stabilized, physical operations. Additionally, methods to minimize the 
number of physical rotations required are detailed. 

As discussed in section [HI by adjusting the phase of our Raman gate laser, we can 
perform rotations about any axis in the equatorial plane of the Bloch sphere. Dehning 
the angle that axis makes with the x-axis as 93 , our physical operations R^{ti/2) are 
composite 7r/2 rotations stabilized with the PBl sequence. While we have constructed 
various analytic solutions to perform arbitrary unitaries in a fixed number of axis angles 
(pj, direct numerical search has proved more efficient. In a test of 1000 randomly 
generated unitaries, we always hnd a solution in three or four rotations, with an average 
of 3.25. 

Our compiler takes a set of requested circuit-level gates for the ion chain and 
computes a target unitary Ui for each ion i. First, Uq is performed on the 01 ion 
pair, using the rotations calculated from our numerical search. While this leaves ion 0 
in the desired state, the physical operations act in a 2 ;-rotated basis due to the phase 
front curvature of the Raman gate beams (figure [2]). Each physical rotation on ion 1 is 
also through an angle of 7r/2; the rotation due to the crosstalk on ion 2 is much smaller. 
Neglecting crosstalk beyond next-nearest neighbors, the resulting physical operations 
are summarized in step 1 of table lAli where the net physical unitary denotes the 
physical effect on ion j of a unitary targeting ion i. In each step, the unitary is 

performed, while the compiler stores a unitary corresponding to the calibrated crosstalk 
on j = [i P IQ + 2]. In step 2, we then need to perform the product of the targeted 
unitary Ui and the inverse of the stored unitary from the previous step. Thus, in each 
subsequent step, we numerically search for physical rotations whose net effect is to undo 
the crosstalk from previous operations and perform the targeted unitary f/j. 

To correct crosstalk, we store a unitary corresponding to the history for each ion 
throughout the algorithm. However, we can take further advantage of this storage to 
reduce the number of applied physical i?(^(7r/2), thereby improving overall hdelity. This 
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Table Al. In contrast to the idealized pairwise control shown in figure [TJb), here we 
show how our compiler must account for crosstalk. 


Addressed 




Step Pair 

Ion 0 Ion 1 

Ion 2 

Ion 3 

1 01 

= !/„ uP 

uP 

- 

2 12 

up = U, 

{CfuP 


3 23 

~ ~ 

u? = U2 


4 34 

- - 

- 

= Us 

Net gate: 

Uo U, 

U2 

Us 


is implemented in two ways. First, subsequent single-qubit gates are concatenated for 
each ion, whenever possible. If they do not need to be performed, we append the inverse 
operation to the stored unitaries. At the next physical operation cascade, the inverse 
history is composed with any additional requested gate, allowing us to perform multiple 
single qubit gates and crosstalk corrections in a single cascade of operations. Second, 
we often perform requested unitaries up to a stored Rz- An arbitrary unitary can 
almost always be performed in only two R^{tt/2) rotations, up to a missing Rz, which 
is nearly half that required for an exact unitary. This has the effect of transforming 
an ion’s relative frame and can simply be appended to the stored unitary. Subsequent 
operations on this ion need only shift the Raman gate phase to account for the frame 
shift. In fact, without global operations, there is no need for a common dehnition of 
zero phase. Instead we dehne a phase zero for each ion by the Raman beam phase front 
when that ion is targeted. 

In some instances, performing unitaries up to a phase still produces an exact result. 
For example, state preparation on the hducial state with every ion in |0) is the same 
with an Rz inserted: 

t/e.act |0) = (l7e.actRi) Rz |0) = U'Rz |0) = U' |0) (A.l) 

since z-rotations do not change |0). Similarly, operations immediately prior to detection 
can be performed up to an Rz without effecting the probability to be in |0) or |1). 

Prior to performing an MS gate, the stored unitary of all targeted ions must be 
physically undone as arbitrary single-qubit operations do not generally commute with 
the MS gate. Additionally, the relative frames of targeted ions must be rotated to match 
the MS gate beams. It would seem we are forced to perform exact single-qubit unitary 
operations around all MS gates. However, these requirements can be softened. If we 
diagonalize the MS interaction with appropriately rotated Hadamard gates on each ion, 
Rz operations will commute and can be tracked to the other side. This is the origin of 
the “gate optimization” additions to hgure EJ^a) and results in only two R^in/2) prior 
to the MS gate and three after it. Exact unitaries would require six to eight R^{'k/2) 
operations. 
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In summary, our compiler stores a unitary history for operations on each ion. 
This enables us to track crosstalk corrections, unperformed single-qubit gates, and 
unperformed Rz rotations. We are able to program algorithms directly from a quantum 
circuit, and the compiler produces a compact set of operations which limit the number 
of physical rotations actually performed. 

Appendix B. MS Echo 

We describe an “echoing” scheme that enables pairwise entanglement despite the 
imperfect pairwise addressing of our MS gate beams. The intuitive motivation is that 
flipping the state of an ion in the MS basis (i.e. -f-)- |—) 3 ,) halfway through the gate 

will cause the net effect of the MS interaction on this ion to be identity. In fact, we 
show that the scheme is exact (within the approximations outlined below) and could 
be used to entangle any m ions in an N ion chain with single-qubit flips on the 

untargeted ions; these “echo” operations must be equally distributed throughout the 
MS gate, requiring it to be performed in equal steps. 

Following |30], we write the interaction picture Hamiltonian, ifint, for two Raman 
beams with symmetric detuning 6 from the carrier transition, assuming identical 
intensity prohles for the red- and blue-detuned beams. Assuming the detuning is close 
to one sideband frequency z/, we ignore all other modes as well as terms oscillating at 
V + 5. In the Lamb-Dicke approximation, is given by 

Hint ~ \/2VLr]Jx [x cos(z/ — 5)t -|- psin(z/ — S)t ], (B.l) 

where H and rj are the maximum Rabi rate and Lamb-Dicke factor, respectively, for any 
ion in the chain, x is the position operator and p the momentum operator (dehned in 
[30]). t is time, and we have dehned 

^ IQ.7?. 

Jx = y^ CiXi and q = - 77 —, (B. 2 ) 

where Xj is the Pauli X operator acting on the ith ion. The Rabi rate Dj and Lamb- 
Dicke factor rji are different for each ion i in the X-ion chain. Equation fIB.ll) also 
neglects off-resonant carrier transitions, which is a valid approximation provided Q 6 
[30] . This Hamiltonian is a generalization of equation ( 6 ) in [30], and the results for the 
propagator U follow directly with a modihed spin-dependence, Jx- If all ions had the 
same Rabi rate and Lamb-Dicke factor, then all q would be | and Jx would reduce to 
the total angular momentum operator Jx- 

The authors of |30| describe how the propagator resulting from a Hamiltonian of 
the form fIB.ip causes eigenstates of Jx to traverse a circular path in phase space. For 
a propagation time 27m/(z/ — 6) with n = 1,2,3..., the phase space trajectory returns 
to the origin, and the ions are left in the original motional state. At these times, the 
propagator U = exp(—iAJ,^), where A is proportional to the geometric area enclosed by 
the loop in phase space. In our general case, the proportionality factor is now a sum 
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over pairs of X operators, weighted by the relative intensity and Lamb-Dicke factor [cj, 
dehned in flB.2|) ]. Simplifying, 

_ N N N N N 

= CiC.XiXj = E c? + 2 E E CiCjX.Xj (B.3) 

2 = 1 j = l 2=1 2 = 1 j>i 

where the first term on the right-hand side is a global phase, independent of the internal 
states of the ions; we ignore this global phase and redefine 


N N 

t/ = exp ( -2iA Y CiCjXiXj ) . 

2=1 j>i 


We define an “echo” operation 




Ek = exp ( ) = -ihfc 


(B.4) 


(B.5) 


which swaps x-eigenstates of the fcth ion through a tt rotation about the y-axis. Unitary 
transformations of the propagator are built from the 


U, 


kik2---kL 




Kl=l 


Kl=l 


= n -2i2l E E ‘:XiX,Xj n 25 . 




2=1 j>i 


h j (B.6) 

JV JV \ /_ \ 

^ (B.7) 

= exp (B.8) 

(B.9) 


where flB.811 can be shown by expanding the operator exponential with a Taylor series, 
and flB.gp uses XjUfe = {—lY^’^Y^Xi. 

We assert that a subspace D within an iV-ion chain can be decoupled so that only 
the ions not in D participate in the net MS interaction through products of these Uk^... 
with ki G D. For example, the following sequence will decouple the kth ion in a Wion 
chain 


N N 

nd,)vA'q 

( \ 

-2iA Y CiCj 1 

1125*, 

2=1 j>i 

\ i / 

V 1 / 

( N N 

- 

1 

1 -2141 


X,X, 

[ i=l j>i 

. 1 

J 


N N 

UkU = exp <1 -2iy4 Y CCj 

2=1 j>i 


l + (-l) 




X,Xj 


(B.IO) 


Since j ^ i, k will will only equal i or j or neither. Therefore, if i or j is equal to fc, the 
term in (IB.lOp in square braces vanishes; conversely, the terms where k ^ i,j survive, 
exactly as desired for an MS gate on all ions but the fcth. 

Generally, a sequence to decouple an d-qubit subspace D can be constructed from 
the product of the Ukik 2 --- where the ion labels ki are all combinations of 0 to d indices 
in D. Table IBTI summarizes possible decoupling sequences up to d = 4. The number of 
such operators is given by a sum over binomial coefficients: 
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Table Bl. Sequences to decouple up to d = 4 ions. The transformed U operators 
all commute. In practice, there may be an optimal order in which to perform these 
operations. 

d Possible sequence 

1 UkU 

2 UuUkUiU 

3 UklmUklUkmUlmUkUlUjnU 

4 Ukl mn Ukl m^kln Ukmn ^Imn UklUkmUknUlmUInUkUI U^U^U 


For an MS gate on m ions out of an N ion chain there are d = N — m ions in D, requiring 
the MS gate be performed in equal steps as initially asserted. Finally, we note that 
since the transformed U operators all commute, they can always be rearranged to give 
a decoupling sequence with a single Ek between each MS interaction. Generally, this 
implies the ability to optimize the echo sequence for any single-qubit control scheme. 

In section 14.21 we demonstrate the most rudimentary echoing scheme. By 
decoupling one ion within a three-ion chain, we generate pairwise entanglement despite 
imperfect pairwise addressing. Within a longer chain, we note one could shorten the 
required echo sequence by neglecting to echo ions which have a near-zero interaction 
(cj) with the gate beams. 

While this echoing technique can decouple any arbitrarily connected subset of 
ions, we have only considered a single vibrational mode under the assumption of near¬ 
detuning. Spectator modes, which tend to lead to residual motional entanglement, may 
be a key error source. In the future, we anticipate using modulated sequences, like those 
described in [33], to alleviate this issue. Such sequences could also be interleaved with 
echo pulses to arbitrarily decouple selected ions with high hdelity. 
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